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Abstract 

With the availabihty of high frequency financial data, nonparametric estimation of volatihty 
of an asset return process becomes feasible. A major problem is how to estimate the volatility 
consistently and efficiently, when the observed asset returns contain error or noise, for example, 
in the form of microstructure noise. The former (consistency ) has been addressed in the recent 
literat ure. However, the resulting estimator is not efficient. In lZhang. Mvkland. and Ait-Sahalial 
()2005|) . the best estimator converges to the true volatility only at the rate of n~^^^. In this 
paper, we propose an estimator, the Multi-scale Realized Volatility (MSRV), which converges to 
the true volatility at the rate of which is the best attainable. We have shown a central 

limit theorem for the MSRV estimator, which permits setting intervals for the true integrated 
volatility on the basis of MSRV. 
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1 Introduction 



This paper is about how to estimate volatihty non-parametrically and efficiently. 



With the availability of high frequency financial data, nonparametric estimation of volatility 
of an asset return process becomes feasible. A major problem is how to estimate the volatility 
consistently and efficiently, when the observed asset returns are noisy. The former (consistency) 
has been addressed in the recent l iterat ure. However, the resulting estimator is not efficient. In 
Zhang. Mvkland. and Ai't-Sahalia (|2005l ^ . the best estimator converges to the true volatility only 
at the rate of n~^/^. In this paper, we propose an estimator which converges to the true volatility 
at the rate of which is the best attainable. The new estimator remains consistent when the 

observation noise is dependent. We call the estimator the Multi Scale Realized Volatility (MSRV) 

To demonstrate the idea, consider {Y} as the observed log prices of a financial instrument, and 
the observations take place at the grid of time points Qn = {tn,i,i = 0, 1, 2, • • • n} that span the time 
interval [0, T]. For the purposes of asymptotics, we shall let Qn become dense in [0, T] as n ^ oo. 

Suppose that {It^^} are noisy, the corresponding the true (latent) log prices are {X}. Their 
relation can be modeled as, 

Yt - = Xt .+et .. (1) 
where t^.i G Gn- The noise ej^ -s will be assumed to be independent of X and iid. 

The mode l in dH) i s quit e realistic, as e videnced by the existence of microstructure noi se in the 
price process (|Brownl (|l99nh . Izhoul (|l99d ^. Icorsi. Zumbach. Muller. and Dacorognal (|200ll )^. 



We further assume that the true log prices {X} satisfy the following equation: 



dXt = 

where Bt is a standard Brownian motion, 
coefficient at are stochastic in the sense that 



: fitdt + atdBt (2) 
Typically, the drift coefficient fit and the diffusion 



dXt{uj) = fi{t,uj)dt + a{t,uj)dBt{Ld) 



(3) 



Throughout this paper, we use the notation in Q to denote (jSJ- By the model in ©, we mean 
that {X} follows an Ito process. A special case is that {X} is Markov, where fit = fj.{t,Xt), and 
at = a{t,Xt). In financial literature, at is called the instantaneous volatility of X. 

Our goal is to estimate Jq atdt, where T can be a day, a month, or other time horizon(s). For 
simplicity, we call afdt the integrated volatility, and denote it by 

{X,X) = r a^dt. 
Jo 
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The general question is, how to estimate nonpar ametrically crfdt, if one can only observe the 
noisy data Yf^. at discrete times G ^n- Gn is formally defined in Sectional 

To the best of our knowledge, there are two types of nonparametric estimators for a^dt in 
the current literature. The first type, the simpler one, is to sum up all the squared returns in [0, T]: 



[Y,Y] 



Y, 



(4) 



this estimator is generally called realized volatility or realized variance (or RV for short). However, 
it has b een reported t hat realized volatility using high-frequency data is not desi rable (see, for 
example. Brown ( 199Cll ) , Zhou ( IQQsh , Corsi. Zumbach. Muller. and Dacorognal ( 2001 ) ) . The reason 
is that it is not consistent, even if the noisy observations Y are available continuously. Under discrete 
observations, the bias and the variance of the realized volatility are the same order as the sample 



size n. 



A slight modification of (jlj is to use the sum of squared returns from a "sparsely selected" 
sample, that is, using a subgrid of Gn- The idea is that by using sparse data, one reduces the 
bias and the variance of the conventional realized volatility. This approach is quite popular in the 
empirical finance literature. However, this "sparse" estimator is still not consistent, in addition, 
which data to subsample and which to discard is arbitra ry. The behavior of this type of est i mator , 
and a sufficiency based improvement of it, is analyzed in Izhang. Mvkland. and Ai't-Sahahal (|2005h . 



A second type of estimator for a}dt is based on t wo sampling scales. As introduced in Section 
4 (p. 1402) of lzhaug. Mvk1a,nd. a,nd AVt-Sa,ha,lial the Two Scales Reahzed Volatility (TSRV) 

has the form 



(TSRV) I n- K +1 / 1 N 

{x,x) = [y,y]("'^)-2 " [y,y]("-i). 



where 



[y,y] 



(n,K) 



1 



nK 



(5) 



(6) 



with K being a positive integer. Thus the estimator in ^ averages the squared returns from sam- 
pling every data point ([y, y]^'^^) and those from sampling every K-ih data point {\Y,Y]^'^^). 
Its asymptotic behavior was derived when K — > oo as n oo. The TSRV estimator has 
many desirable features, including asymptotic unbiasedness, consistency, and asymptotic normal- 
ity^. However, its rate of convergence i s not satisfactory. For an instance, the best estimator in 



Zhang. Mvkland. and Ai't-Sahalial (|200fJ ) converges to Jq a^dt at the rate of 



n 



-1/6 



In this paper, we propose a new class of estimators, collectively referred to as Multi Scale 
Realized Volatility (MSRV) which converge to afdt at the rate of n~^/^. This new estimator has 

related estimator can be found in lZhouTilQQSl) and lHansen and Lundel i2006l) . however, their estimator (takes 
k to be fixed) does not yield a consistent estimator. 
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the form, 



in) 



M 

E 

1=1 



ai[Y,Y] 



{n,K,) 



(TSRV) 

where M is a positive integer greater than 2. Comparing to {X,X)rp which uses two time 

(") 

scales (1 and K), {X,X) combines M different time scales. The weights Oj are selected so that 

(n) 

{X, X) is unbiased and has optimal convergence rate. The rationale is that by combining more 
than two time scales, we can improve the efficiency of the estimator. Interestingly, the n~^/^ rate of 
convergence in our new estimato r is the same as the one in parametric estimation for volatility, when 
the true process is Markov (see Gloter and .TacodI (|2nnril ) ) . Thus this rate is the bes t attainable . 
Earlier related results in the same direction can be found in IsteinI Jl98ill99dll99^ and lYinal Jl99l] . 
199,lh . S ee also |A it-Sahalia. Mv kland. and Z hang* ('2005a[L 

Relate d independent work can also be 



found in 



Barndorff-Nielsen. Hansen. Lunde. and Sheohardl ll2004 ). For the est imating functions- 



based approach, there is a nice review by Bibbv. Jacobsen. and S0rensen ( 2002^ . 



We emphasize that our MSRV estimator is nonpar ametric, and the true process follows a more 
general Ito process, where the volatility could depend on the entire history of the X process plus 
additional randomness. 

The paper is organized as following. In section 2, we motivate the idea of averaging over M 
different time scales. As we shall see, our estimator is unbiased, and its asymptotic variance comes 
from the noise (the Cf^jS) as well as from the discreteness of the sampling times In Sections 
3-4, we derive the weights Oj's which are optimal for minimizing the variance that comes from noise, 
and we give a central limit theorem for the contribution of the noise term. A specific family of 
weights is introduced in section 4. We then elaborate on the discretization error in Section 5, and 
show a CLT for this error. Section 6 the gives the central limit theorem for the MSRV estimator. 

For the statements of results, we shall use the following assumptions: 

Assumption 1. (Structure of the latent process) . The X process is adapted to a filtration (^Yj), and 
satisfies where Bt is an {Xt)-Brownian motion, and the /Uj and at are {Xt)-adapted processes 
which are continuous almost surely. Also both processes are bounded above by a constant, and at is 
bounded away from zero. We denote X = Xt. 

As a technical matter, we suppose that there is a a-field J\f and a continuous finite dimensional 
local martingale (Mt) so that Xt = a{Ms, < s < t)\/ J\f . 

Assumption 2. (Structure of the noise). The e^^ . are independent and identically distributed, 
with E[e]=Q and E[e^] < oo. The et„ j are also independent of X 



These assumptions a re not minimal for all resu lts. In terms of the structure of the proce s s, see , 



for example. Section 5 in lJacod and Protteil ()l998l ^ and Proposition 1 in iMvkland and Zhand \2m'h 

for examples of statements where the (j, and a processes are not assumed to be continuous. For the 
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method ology to incorporate dependence into the noise structure, see lAi't-SahaUa. Mvkland. and Zhane 
( 2nn5bl 'l. Our current assumptions, however, provide a setup with substantial generality without 
overly complicating the proofs. 

The final item in Assumption [T] is s t andar d for the type of limit result that we discus s, cf. 
similar conditions in 



.Tacod and 



Zhang. Mvkland. and Ai't-Sahalia 



rotteH (|199HD . IZhand (|2nnih . iMvkland and Zhand (j2QQa) and 



2 Motivation: Averaging the Observations of {X, X) 



In 



Zhang. Mvkland. and Aiit-Sahahal l|2nn,4 ). we have observed that by combining the square incre- 

~ (TSRV) 

ments of the returns from two time scales, the resulting two-scale estimator (X, X)rp in ^ 

improves upon the realized volatility, which uses only one time scale, as in @. The improvement 
is about reducing both the bias and the variance. 

If the two-scale estimator is better than the one-scale estimator, a natural question would be 
how about the estimator combining more than 2 time scales. This question motivates the present 
paper. In this section we briefly go through the main argument. 



To proceed, recall definition © of [y, y]*^"'^), and set, similarly, 



{n,K) 



1 



and 



tn,i&Qn,i>K 

Under one can decompose [y, y]^"'''^'* into 

[y,y]("'^) = + [e,e]('"-^) + 2[X,e]("'^). 



(7) 



We consider estimators on the form 



(n) 



M 



{x,x) = j;ai[y,y]("^ 



i=l 



(8) 



where a^'s are the weights to be determined. A first intuitive requirement is obtained by noting 
that 

E{{Xxy)x process ) = ^ ai[X, + 2Ee' ^ + ' 

i=l 4=1 



(9) 



Multi Scale Realized Volatility 



5 



Since [X, are asymptotically unbiased for {X, X) (|zhang. Mvkland. and Ait-Sahalial (jiooi)), 

it is natural to require that 

(10) 



Oi = 1 and ai — ^—rz = 



1=1 1=1 
A slight redefinition will now make the problem more transparent. Let 



ai = «! 



(n + 1) 



1 1 



'1 



Our conditions on the a's are now equivalent to 
Condition 1. = 1, 



02 = 02 — (ai — ai) and a-i = for i > 3. (11) 



Condition 2. ^fi^ = 



(n) 



To understand the estimator {X,X) in terms of the aj's, consider the following asymptotic 
statement. Here, and everywhere below, we allow Oj, Ki and M to depend on n (i.e., they have 
the form a„ j, j and M„), though sometimes the dependence on n is suppressed in the notation. 
We obtain (for proof, see Section [HJ 



Proposition 1. Suppose that K^^i and Kn,2 o,Te 0(1) as n ^ oo. Under Assumptions^S^ 

M 



^ ^ (n) 

{X,X) =^a,[y,y]("'^') -2Ee2 + Op(n-i 



-l/2> 



(12) 



1=1 



To further analyze the terms in ()12() . write 



[Y, = [X, + ^Y1 + ^-'^ + ^-.^ 

where Un,K will turn out to be the main error term, 

n 



(13) 



_ 2 



K 



(14) 



and Vn,K will be a remainder term, given by Vn,K = 2[X, ej^"'-^) - Ej^o^ ~ :^ EUn-i^+i 
We now can see the impact of Condition [21 To wit, from equation H12|) . 



(n) 



M 



M n 

a. 



M 



M 



{X,X)" = ^a,[X,X]("'^>) + 2^-i^6?„,^+^a,[/„,x.+E«*^-.^.-2^^' + Op(^"' 



i=l j=0 



i=l 



M 



M 



n 



-1/2N 



(15) 
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where Rn is the overah remainder term, i2„ = "^fLi CLiVn,Ki — 2Ee^. Thus, apart from the con- 
tribution of this remainder term, Condition [21 removes the bias term due to X^e^j, not only in 
expectation, but almost surely. We emphasize this to stress that though we have assumed that 
the et^ . are i.i.d., our estimator is quite robust to the nature of the noise. As before, Condition 1 
assures that the first term in H15() will be asymptotically unbiased for {X, X) . 

Furthermore, for i / /, the Un,Ki and Un,Ki are uncorrelated. Since Un^Ki and Un^Ki are also 
the end points of zero-mean martingales, they are asymptotically independent as n ^ oo. Finally, 
the last term Rn is treated separately in the proof of Theorem ^ For now, we focus on the terms 
other than the Vn,Ki 's. 

If one presupposes Condition [21 and that i?„ is comparatively small, it is as if we observe 

[X,X](^') + C/„,^,, i = l,...,M. 

Under the ideal world of continuous observations (that is, if we take [XjX]^^"-^ to stand in for 
{X,X)), Condition [21 makes it possible that we get M (almost) independent measurements of 
{X,X). This motivates the form of the MSRV estimator. 

Our aim is to use Conditions I1I2I to construct optimal weights Oj. We proceed to investigate 
what happens if we just take [X, X]^^*) ~ {X,X) in Section 3-4. From Section 5 on, we consider 
the more exact calculation that follows from = {X,X) + Op{{n/ Ki)-'^/'^). 



3 Asymptotics for the Noise Term 



As above, to get a meaningful asymptotics, we let all quantities depend on n, thus Oj = a„_j, 
M = Mn, Ki = Kn,i, [Y, Y]^^'^ = [Y, etc. Sometimes the dependence on n is suppressed in 

the notation. All results are proved in Section [HI 

Consider first the noise term 

Cn = X] an,iUn,K„,i (16) 
i=l 

The variance of C„ is as follows. 

Proposition 2. (Variance of the noise term.). Set — ^Yli=i ("^^)^- Suppose that the ef^ - are 
iid, with mean zero and Ee^ < oo, and that Mn = o(n) as n — > cxo. Then 

Var{Cn)=llniEe''f (1 + 0(1)). (17) 

Also, jn minimized, subject to Conditions by choosing 

I^n,i(Kn,i Kn) 



MnVar(Kn 



(18) 
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where Kn = ^ J^fti Kn,i and Var{Kn) = ^ E»=l ^l^i - (i E*="i Kn,if ■ The resulting mini- 
mal value of 'jn is 

*2 4 

M„yar(/^„)- ^^^^ 



Sin ce the Un,K are end points of martingales, by the martingale central limit theorem (|Hall and Hevde 
Chapter 3), we obtain more precisely the following: 



Theorem 1. Suppose that the et„. are iid, with Ee^ < oo, and that M = Mn = o{n) as n ^ oo. 
Suppose that maxi<j<Af„ \an,i/ {i^n)\ —>■ as n ^ oo. Then Cn/ {n^^'^ln) N(0, E^e^)"^) in law, 
both unconditionally and conditionally on X . ■ 



4 A Class of Estimators, and Further Asymptotics for the Noise 
Term 

We here develop a class of weights an,i which we shall use in the rest of the paper. The precise 
form of the weights is given in Theorem |5J The rest of this section is motivation. 

In the following and for the rest of the paper, assume that all scales i = 1, M are used, which 
is to say that Kn^i = i- In this case, Kn = (M„ + l)/2 and Var{Kn) = (M^ — 1)/12, and the 
optimal weights from Proposition |3 are then given by 

ani = 12 ^ ^ " (20) 

The minimum variance is given through 7*^ = 48/[M„(M^ — 1)], so that 

Var{Cn) = ^Sn{Ee^ f /[Mn{Ml - 1)]. 
The form (|2U|) motivates us to consider weights on the form 

an, = ^i.M„(^), i = l,...,M„, (21) 

as this gives rise to a tractable class of estimators. We specifically take: 

wm{x) = xh{x) + Ar'^xhi{x) + M~'^xh2{x) + M-^xh3{x) + o(M~^), (22) 

where h and hi are functions independent of M. The reason for considering this particular func- 
tional form, where WMix) must suitably vanish at zero, is that condition ((2)) translates roughly into 
a requirement that ^^^^j^dx be approximately zero. 

In terms of conditions on the function h, Conditions (EO)-® imply that we have to make the 
following requirements on h: 
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Condition 3. xh{x)dx = 1, 
Condition 4. h{x)dx = 0. 

With slightly stronger requirement on h, we can show that (|15|) holds more generally. 

Theorem 2. Let ho = h, and suppose that for i = 0, 2, hi is 3—i times continuously differentiahle 
on [0,1], and that /13 is continuous on [0,1]. Suppose that h satisfies Conditions\^l^ Also assume 
that 

"1 1 

hi{x)dx + -{h{l) - h{{))) = 0, 

j\^{x)dx + ]^{hr{l)-him + ^{h'{l)-h'm = 0, (23) 

1 

and / h-i{x)dx + —{h[{l) - h[{{))) = 0. 

Let the a„^j 6e given by \21\) - 1i^) . where the o(M~'^) is uniform in x £ [0,1]. Finally, suppose that 
the et^- are i.i.d., with Ee^ < 00. Then approximation kL5]) remains valid, up to 0p{n/M'^). 

The final class of estimators. Our estimation procedure will in the following be using weights 
o„^j which satisfy the description in Theorem [2 

Remark 1. [Comments on Theorem\^] By adding terms in \2^) . one can make the approxima- 
tion in \15\) as good as one wants (up to Op{n~^^'^)). We will later use Mn = 0{n^^'^), which is 
why we have chosen the given number of terms in h2'2\) . Also, it should be noted that the approx- 
imation to Condition has to be much finer than to Condition Q since we are seeking to make 
YhLi jti TJi=o = n (riiii ^) + Op(l)) negligible for asymptotic purposes. 

As we shall see, the specific choices for hi, /i2, and /13 do not play any role in any of the later 
expressions for asymptotic variance. A simple choice of hi which satisfies (j23|) is given by hi{x) = 
—h'{x)/2, with h2{x) = /12 and h^{x) = h^, both constants. In this case, /i2 = — (/i'(l) — /i'(0))/6 
and hs = (/i"(l) - /i"(0))/24. With this choice, one obtains 

For the noise-optimal weights in (|2()j) at the end of Section 3, h takes the form 

hl{x) = 12 (^x-^y (25) 

Under this choice, the an,i given by PH) is identical to the one in (^0]) . up to a negligible multi- 
plicative factor of (1 - M~^)~^. 



The following corollary to Theorem ^ is now immediate, since 7^ = 4M-^ £h{xfdx{l + oil)) 
as n — > 00. 
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Corollary 1. Suppose that the ej^ - are iid, with Ee^ < oo, and that M = Mn = o{n) as n ^ oo. 
Also assume that the On^i are given by \21\) . and that the conditions of Theorem are satisfied. 
Then {M^/n)^^'^(n N{0,4:E{e'^)'^ h{x)'^dx in law, both unconditionally and conditionally on 
X. I 



5 Asymptotics of the Discretization Error 

We have obtained the optimal weights as far as reducing the noise is concerned. However, as in 
H15|l . there remains two types of error: the discretization error., due to the fact that the observations 
only take place at discrete time points, along with the residual Rn, which also will turn out to not 
quite vanish. We study these in turn, and then state a result for the total asymptotics for the 
MSRV estimator. 

For the discretization error, we need some additional concepts. 

Definition 1. Let = tn,o < tn,i < ■■■ < tn,n = T be the observation times when there are 
n observations. We refer to Qn = jtn.n. i , .... tn. w) as a, "grid" or a "partition" of [0, T]. Fol- 



lowing Section 2.6 of Mvkland and Zhanci 1 20 Oi ). the "Asymptotic Quadratic Variation of Time" 



("AQVT") Hit) is defined by 

H{t) =lim^ V {tn,i-tn,i-lf, (26) 



provided the limit exists. 
We assume that 

max \tn,i+i - tn,i\ = O ( - ) , (27) 
l<i<n I ' ^ ' I \n J 

whence every subsequence has a subsequence so that the asymptotic quadratic variation of time 
exists. From an applied point of view, there is little loss in assum ing the existence of the asvmptotic 
quadr atic variation of time, cf. the argument at the very end of Zhang. Mvkland. and A'it-Sahalial 
(|2nn,4 l (on p. 1411). 

Note that from (|27() . H{t) is Lipschitz continuous provided it exists. We give the following 
change-of-variable rule for the AQVT: 



Lemma 1. (Change of variables in the AQVT.) Assume \2'T\) and that the AQVT H(t) exists. Let 
G : [0, T] [0, T] be Lipschitz continuous. Set Un,i = G{tn,i)- Then 



K{u) = lim — (n„ i - n„ 

n— >oo 1 '—^ 

U„A<U 
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exists, and 

H'{t)G'{t) = K'{G{t)) (28) 

almost everywhere on [0, T] . 

The following result is also useful and illustrative. 
Lemma 2. Assume the conditions of Lemma^ Then K(T) = T if and only if 

( Un,i - Un,i-i ) = o(n"^). (29) 

Remark 2. The importance of these two lemmas is that one can compare irregular and "almost 
equidistant" sampling. If H'{t) exists, is continuous, and is hounded below by a constant c > 0, 
one can define G{t) = Jq H'{s)~^ds, and consider the process Xu = Xg{u)- This process satisfies 
the same regularity conditions as those that we impose on X, and, furthermore, the sampling times 
Un.i = G{tn,i) are close to equidistant in the sense of equation h2y\) . The further implication of this 
is discussed in Remark after Theorem 

Define ry as the nonnegative square root of 

rp 

rf=f H'{t)afdt (30) 

Finally, we define "stable convergence". 







Definition 2. If Zn is a sequence of -measurable random variables, the Z„ converges stably in law 
to Z as n ^ oo if there is an extension of X so that for all A ^ X and for all bounded continuous 
g, ElAg{Zn) ElAg{Z) as n ^ oo. 



For further discussi on of stable con v ergen ce , see lEenvil Jl963h . lAldous and Eaelesonl (119781) . 



Chapter 3 (p. 56) of iHall and He^ Jl98nl 'l. lEootzenl (|l98nl ) and Section 2 (p. 169-170) of 



Jacod and Protten ()l998l ). It is a useful device in operationalizing asymptotic conditionality. There 



is some choice in what one takes as the u-field X in this definition. 

We can now state the main theorem for the asymptotic behavior of finitely many of the 

Theorem 3. (CLT for the discretization error in [X,X\^^\) Suppose the structure of X follows 
Assumption^ Also suppose that the observation times tn^i are nonrandom, satisfy \2T\j , and that 
the asymptotic quadratic variation of time II{t) exists and is continuously differentiable. Assume 
that mino<t<T -f^'(i) > 0. Let Mn oo as n ^ oo, with Mn = o(n). Let (Kn^i, Kn^L)/Mn — > 
(ki, hl) as n ^ oo. Let T be an L x L matrix with (I, J) entry given by 

r '^rp . f .fo min(K/,Kj)\ 

r7,j=-rmm(K/,Kj) 3 , (31) 

3 \ max(K7, nj) J 
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and let Z be a normal random vector with covariance matrix T. Let Z be independent of X . Then, 
asn^oo the ?;ector (n/M„)i/2([X, - (X, X) [X, - (X, X)) converges stably 

in law to rjZ . 

Remark 3. Even in the scalar fL = 1) case, this result in Th eorem\^is a gain over our earlier The- 
orem 3 (p. 1401) in Zhang. Mukland. and Ait-Sahalid 1200,^) . To characterize the asymptotic dis- 
tribution we use an asymptotic quadratic variation of time fA QVT) which is independent of choice 
of scale and coincides with the original object introduced in \Mvk land and Zhanc^ \200i ) (Section 
2.6). This is unlike the time variation measure used in section 3.4 in Zha na. Mvkland. and Ait-Sahali^ 
\200d) . and Theorem provides a substantial simplification of the asymptotic expressions. To do 
this, we have used the approach described above in Remark\^ 



It is co njectured that the repularitu conditions for Theorem can be reduced to those of Propo- 
sition 1 of Mykland and Zho,nc \200^) . but investigating this is beyond the scope of this paper. 



As a corollary to Theorem |21 we now finally obtain the asymptotics for the discretization part 
of the MSRV, as follows. 

Corollary 2. (CLT for the discretization error in the MSRV.) Let an^i satisfy \21\) - \2'J\) . and let 
the conditions of Theorem be satisfied. Further, make Assumption Also suppose that the 
observation times tn,i are nonrandom, satisfy \21\j , and that the asymptotic quadratic variation of 
time H{t) exists and is continuously differentiable. Assume that m.m.Q<_t<T H' [t) > 0. Let Mn — > 00 
as n ^ 00, with Mn/n = o(l) and M^/n — > 00. Set 

Vh=o'^V^ Hy)Hx)y'^ {^x - y) dy (32) 

•J Jo Jo 

Then 

(n/M„)V2 f ^a„,i[X,X]("'^) - {X,X)j ^ %Z (33) 
stably in law, where Z is standard normal and independent of X . 

Remark 4. Note that the condition M'^/n 00 is present because we have not imposed too many 
conditions on h; if it were necessary, the assumption could be removed by considering a slightly 
smaller class of hs. 



6 Overall Asymptotics for the MSRV Estimator 



There are two main sources of error in the MSRV. On the one hand, we have seen in Corollary ^ (at 

~ {«) 

the end of Section Q that if Af„ time scales are used, the part of {X,X) — {X,X) which is due 
purely to the noise e can be reduced to have order Op(n^/^Af„ ^^^). At the same time. Corollary 
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l^lshows that the pure discretization error is of order Op{n ^/"^Mn"^). To balance these two terms, 
the optimal M„ is therefore of the order 

M„ = 0(ni/2), (34) 

assuming that the remainder term in H15|) does not cause problems, which is indeed the case. This 
leads to a variance-variance tradeoff, and the rate of convergence for the MSRV estimator is then 
{X,X) - {X,X) = Op(n-V4). This result is an improvement on the two scales estimator, for 
which the corresponding rate is Op(n^^/^). We embody this in the following result. 

Theorem 4. Let an,i satisfy h21]) - Ji^) . and let the conditions of Theorem\^be satisfied. Further, 
make AssumptionsUS^ Also suppose that the observation times tn^i are nonrandom, satisfy i27[ ), 
and that the asymptotic quadratic variation of time H[t) exists and is continuously differentiahle. 
Assume that mino<t<T -f^'(^) > 0. Suppose that Mn/n^l"^ c as n ^ oo. Let Z he a standard 
normal random variable independent of X . Set 

vl = Ac-^{Ee^y j h{xydx + c^Tr]'^ I dxl h{y)h{x)y^ {3x - y) dy 
] Jo Jo 

^1 py pi pi 



4c-^{Ee^f [ h{xfdx + c\Trf [ , 
Jo 3 Jq 

+ 4c-Var(e^) / / xh{x)h{y)dxdy + Scr^Ee"^ 

Jo Jo Jo Jo 



h{x)h{y)mm{x,y)dxdy {X, X) (35) 



Then 

(^(X;^)^"^ - ^ u^,Z, (36) 

stably in law, as n ^ oo. 

For the noise optimal /i- function from equation H25|) (cf. equation ()2U|)'). we can now calculate 
the value of the asymptotic variance of the MSRV. Note that if h{x) = 12(x — 1/2), we obtain 

dx I h{y)h{x)y {3x-y)dy = —, 



Jo 



^ f-y 3 

xh{x)h{y)dxdy 



JO 



5' 

1 ri 6 
h(x)h(y)m.m(x,y)dxdy = -. 

5 



'0 JO 

Hence, in this case, the asymptotic variance becomes 



i^l = 48c-HEe'f + ^^Trf + ^^c-^Var{e^) + ^c-^Ee-" {X,X) (37) 
35 5 5 



7 Conclusion 



In this paper, we have introduced the Multi Scale Realized Volatility (MSRV) and shown a central 
limit theorem (Theorem ^ for this estimator. This permits the setting of intervals for the true 
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integrated volatility on the basis of the MSRV. As a consequence of our result, it is clear that the 
MSRV is rate efficient, with a rate of convergence of Op(n^^/'^). 

In terms of the general study of realized volatilities. Section [3 also shows further properties of 



the as ymp totic quadratic variation of time fAQVT) , as earlier introduced bv lMvkland and Zhang 



and 



Zhang. Mvkland. and Ai't-Sahalial pOOfj ) . In particular, Theorem 01 shows that one can 



use the regular one-s tep AQVT also for multistep realized vol atilities, thus improving on Theorems 



2 and 3 (p. 1401) in lZhang. Mvkland. and Ai't-Sahalia 



Finally, note that most of the arguments we have used hold up also when the noise process et^ . 
is no longer iid. One can, for example, model this process as being stationary (but with mean zero). 
If the process is sufficiently mixing, this will change the asymptotic variance of the MSRV, but not 



the co nsistency, nor the convergence rate of Op{n ^/^), see for example Chapter 5 of Hall and Hevde 
for the basic limit theory for dependent sums. However, we have not sought to develop the 
specific conditions for the CLT to hold in the case when the process is mixing. 

8 Proofs of Results. 

Note that for ease of notation, we sometimes suppress the dependence on n in the notation. For 
example, Oj = a„_j, M = Mn, Ki = Kn,i, [Y,Y]^^') = [Y,Y]^'^'^\ etc. Also, we in this section write 
ti for tn,i, to avoid cluttering of the notations. 

8.1 Proof of Proposition [T], 

Write 

{x,xy = ^a4y,y]('^'^') + (ai -ai)([y,y]('^'^i) - [y,y]("'^2)) 

i=l 
M 

= ^a4y,y]("'^') -2Se2 + Op(n'i/2) (38) 

i=l 

where the final approximation follows from Lemma 1 (p. 1398) in Zhang. Mvkland. and Ai't-Sahalial 
( 200,4 ). 
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8.2 Proof of Proposition m 

Since Un,K„ i ^iid Un,K„ I are uncorrelated (i 7^ I) zero-mean martingales, 

Mn 

Var{Cn) = '^al^iVar{Un,Kr,^,) 



i=l 

M, 



• 1 



n,i 



= jMEe^f (1 + 0(1)), (39) 
showing equation H17|). The last transition in (|39() follows because M„ = o(n). 

We minimize 7^, subject to the constraints in Conditions I1I2I This is established by setting 

g^bl + A,(^ a., - 1) + ^)] =8^ + X, + ^' 



n,i 



to zero, resulting in an,i = —^(XiK^ ^ + A2-f^n,j)- One can determine the A's by solving 
J 1 = On,i = -^(Ai Y^fti Kl i + A2 Y^ti Kn,i) 

\ = E.=l ^ = -|(Ai Ef-il Kn, + MM 

This leads to 

Ai = — — and A2 



where ^„ and yar(i^n) are as given in Proposition This shows the rest of the proposition. 
8.3 Proof of Theorem [H 

Assume without loss of generality that Ki = i for i = 1, M. To avoid cluttering the notation, we 
write ai for an,i- Note that d is the end point of a martingale. We show that Cn/(n^^'^Jn) satisfies 
the conditions of the version of the Mart ingale Central Limit Theorem which is stated in Corollary 
3.1 (p. 58- 59) of Hall an d Hevd3 (|l98ol ). The result then follows. Note that we shall take, in the 



notation of iHall and Hevda (il980|), Tn,j to be the smallest a-field making et^, i = 1, j, and the 
whole Xt process, measurable. 

We start with the Lindeberg condition. For given 5, define fs(x) = E{e'^x'^I^^^^^ygy). Also set 

r„(x) = -— -^et, for ^— <x<L 

V i=i ^ J 



n n 
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We then obtain 



M„Aj 



1 2a,; 



1 " / 1 2a- 

~2_.-^/5ni/2 V ^etj_, 

i=i \ i=i 

/ rn{x)dx since the ej. are i.i.d. 
Jo 



'0 

as n ^ oo, 



(40) 



where the last transition is explained in the nex t paragraph. By Che bychev's inequality, the 
conditional Lindeberg condition in Corollary 3.1 of lHall and Hevdel Jl98nl ) is thus satisfied. 

The last transition in (|4()|1 is because of the following. First fix x G [0,1), and let j„ be 
the corresponding j in the definition of r„(x). Let Z„ = —-^^^i''^" ^^U^ so that r„(x) = 

Note that Z„ is a sum of independent random yariables which, satisfies the Lindeberg condition: 

2 



M„Aj„ 

1=1 



i=l 







as n — > oo, since max, \ ai/i^n\ 0. The ensuing asymptotic normality of Zn (if necessary by going 
to subsequences of subsequences) shows that r^ix) — > as n ^ oo. Since < r„(x) < 1, the final 
transition in (|4()|) follows by dominated convergence. 



We now turn to the sum of conditional variances in the corollary in lHall and Hevdd (jlQSOl). 



1 '^'2a. 



E — 1^,. 



2a,: 



j=i \ i=i 



l + Op(l). 



(41) 



The last transition is obvious by appealing to M-dependence. A rigorous but tedious proof is 
obtained by splitting the sum into main terms of the type e^. and cross-terms of the form et-e^. 

In view of (|4U() - H41|) . T heorem U is proved by using Corollary 3.1 and the Remarks following 
this corollary (p. 58-59) in lHall and Hevdj JlQSoI 'l. 



Multi Scale Realized Volatility 



16 



8.4 Proof of Theorem H 

We need to show that J2i=i ^7=0 = Op{n/M^), in other words, we need Yld=[^ TT^ = 
o(M~^). By Taylor expansion 

1=1 •^'J i=l i=l i=l 

„1 . M . M . 

= / + — (/i(l)-/i(0)) + — — ) — ) + o(M-3) 

-"^ i=i i=i 

= f h{x)dx + ^(/i(l) - /i(O)) + -^i.h'{l) - h'm + o{M-% (42) 



where the later line follows by iterating the first line. By similar argument on hi to /13, 

i=l i=l i=l i=l i=l 

= h{x)dx + ^ (^^' hi{x)dx + ^(/i(l) - /i(O))) 
by ((^ . This shows the result. 



8.5 Proof of Lemma [H 



To get the rigorous statement, we proceed as follows. Every subsequence has a further subsequence 
for which K{u) exists, and this K is obviously Lipschitz continuous. We will show that (|28jl hold. 
Since this equation is independent of subsequence, the result will have been proved. 



Let Bt be a standard Brownian motion, and let Bt = BQ(^iy By comparing the asymptotic 
distributions of (T/n)-V^[y-.^^.fg,.-^,,- < >,] and f r/n)-V^[y-.,^^Ji?„.. -K.- < 

B, B >u\, we obtain from Proposition 1 of lMvkland and Zhand (2003) that 



2H'{s){< B,B >',yds 



Git) 



2K'{v){< B,B>'J^dv for ah t G [O.T]. 



Since < B,B >^= 1 and < B,B >'^= G'{s) a.e., equation (^5]) . and hence the lemma, follows. 
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8.6 Proof of Lemma [51 



n 
T 



Set 5n,i = Un,i - Un,i-i - T/n. Then 

Y{Un,i-Un,i-lf = f Ei (i +^n,i)^ 

Since Yli ^n,i = 0, the Lemma follows by letting n oo. 



8.7 Proof of Theorem 

Following Lemmas n and 121 and Remark |21 we can assume without loss of generality that the tn,i 
satisfy (in place of the Un,i) the equation ((^. 

Consider the scalar case (L = 1) first, with = i^n,i = Afn- In the sequel, all prelimiting 
quantities are subscripted by n, and we suppress the n for ease of notation (except when it seems 
necessary). We now refer to Theorems 2 and 3 (p. 1401) in Zhang. Mvkland. and A'it-Sahalial 
( 2005h . Use the notation Atj, hi and r/„ as i n that paper, and let At = Tin. (Note t hat the usage 



of "ry" in this paper is different from that of lZhang. Mvkland. and Ai't-Sahalial ()2fln5l ) . Also define 

4 



hi 



KAt 



(K-l)Ai 

^ (1 - and 7?^ = J2h. 



,alAt. 



Note that if we show that — ^ in probability as n — > cxd, we have shown the scalar version of 
the theorem. This is because we we will then have shown that the conditions of the two Theorems 



m 



Zhang. Mvkland. and Ai't-Sahalial (|2005l ) are satisfied, and that we can calculate the asymptotic 



variances as if tn,i = iT/n. 
To this end, note first that 



^h,aliAt,-At)\ < (a+)4 (^^hj^ (^^{AU-Atf 
= 0{n^/^) X o(n-i/2) = o(l). 



1/2 



(43) 



where the orders follow, respectively, from equation (45) in Zhang. Mvkland. and Ai't-Sahalia ( 2003 ). 
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and equation H29|) in this paper. Then note that 

K-l 



k)alAi\ = \^{a^)'Y.(^ 



3_ 
K 



K-l , 

i=i ^ 

0(1) X o(n-V2) =o(l) 



«=(ir-i)+ 



4(. 



1/2 



(44) 



where , again, the orders fohow, respectively, from equation (45) in Zhang. Mvkland. and Ai't-Sahahal 
and equation in this paper. 



Equations (|l3|) - (|il|) combine to show that r)„ 



in probabihty as n — > oo. 



For the general (L > 1) case, first note that since /^t and at are bounded (Assumption ITl), b y 
Girsanov's Theorem (see, for exa mple. Chapter 3.5 (pp - 190- 201) of lKaratzas and Shrevd (|199lh . 
or Chapter II-3b (pp. 168-170) of Jacod and ShirvaevI (|2003h 'l. we can without loss of generality 
further suppose that ut = ide ntically. This is because of the stability of the convergence, cf. the 
methodology in Rootzen ( 198o[ ). 

Now set 



n-l 



(X,X) 



{K) 



jA(K-l) 



K 



Y^iXt^,, - Xt^) E (i^-r)(X,^_,,,-X, 



j=o 



r=l 



and note that 



fx, XI ("'-f^) 



(X,X)(^) + + Op{K/n) 

{X, X)f^ + {X, X) + Op(n-i/2) + Op(K/n), 



from Proposition 1 in lMvkland and Zhand |(200j). 



Let M"'^ be the continuous martingale for which M ^'^ = (X, X)^^) (n/M„)^/^. The proof of 



nJ 



Theorem 2 in 



rn,K, 



Zhang. Mvkland. and Ai't-Sahalial (|2005l ^ actually establishes t hat the sequence of 



proces ses (M; ' "'0 is C-tight in the sense of Definition VI. 3. 25 (p. 351) of I Jacod and ShirvaevI 
J20n3h. This is because o f Theorem VI.4.13 (p. 358) and Corollary VI.6.30 (p. 385), also in 



Jacod and ShirvaevI ()2003l ) . The same corollary then establishes that asymptotic distribution is as 
described in Theorem 13 provided we can show that 



>T— > ?7 r as n ^ oo. 



(45) 



This is because of Levy's Theorem (see iKaratzas and Shreve J 1991 ), Theorem 3.16, p. 157). Th e 
stable convergence follows as in the proof of Theorem 3 of Zhang. Mvkland. and Ai't-Sahalia ( 20051 ) , 
the conditions for which have already been satisfied. 
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We finally need to show ()45|) . As in the scalar case, we assume ()29() . and the same kind argument 
used in the scalar case carries over to show that we can take ti^n = iT/n for the purposes of our 
calculation. The computation is then tedious but straigh tforward, and carried out similarly to tha t 
for the quadratic variation in the proof of Theorem 2 in Izhang. Mvkland. and Ait-Sahalial (|2nn,4 ). 
Theorem |31 is thus proved. 



8.8 Proof of Corollary El 

First of all, note that since M^/n oo, Yli=[ '^n,i = o{—{n/Mny/'^). In lieu of equation it is 
therefore enough to prove 

Mn 

{n/Mnf^ an,i {[X, - (X, X)) ^ (46) 

i=l 

Also, as in the proof of Theorem |31 our assumptions imply that we can take fit = identically 
without loss of generality. 

Since there are asymptotically infinitely many [X, Xj'^^'^^'s involved in equation we have to 
approximate with a finite number of these. To this end, let 5 > be an arbitrary number (6 < 1). 
Let a = 1 - 5/V2. Let L be an integer sufficiently large that 2a^ ^ < 6"^. For / = 1,...,L, let 
kj = a^~^ , and kq = 0. For i = 1,...,M„, define li^n to be the value 1,1 < I < L for which 
i/Mn G (k7_i,k/]. Then note that, if ||;7|| = {EU^y/^, 

(n/M„)^/'|| Va„,, - || < (n/M„)i/2 V |a„,,| x max ||[X, 

^ — ' \ / ^ — ' l<i<n 
i=l i=l 

(47) 

Now let in be the value M„ which maximizes given n, 

and let In = Ii„,n- 

For the moment, let be an unbounded set of positive integers so that (i n/Mn, In/Mn)npN con- 



verges . Call the limit {ki,K2)- Bv the proof of TheoremlHland of Theorem 2 in lZhang. Mvkland. and Ait-Sahalia 



( 2nnf)h . (n/M„)([X,X] - [X,X]("'-^"))2 is uniformly integrable. By the statement of Theorem 



13 it then follows that, as n ^ oo through N 

in/Mn)E{[X,X]^^''-^ - [X,X]("'^"))2 ^ Erj\r2,2 + ri,i - 2r,,2) 

= Er,'2T^2 (l-^)' 

< Er]'^T5'^ (48) 

by construction. Since every subsequence has a subsequence for which {in/Mn, In/Mn), it follows 
from equation (|17l) that 

M„ 

limsup(n/M„)i/2|| y^a„ 7[x,A:]("'') - [X, || < 5{ErfTf/'^ max \xh{x)\. (49) 

n^oo , ' V / 0<a;<l 

1=1 
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The result of Corollary |21 thus follows by computing the limit of 

(n/M„)V2^a„,, - {X,X)) , (50) 

1=1 

and then letting S ^ 0. 

8.9 Proof of Theorem m 

The remainder term Rn from equation (|15j) can be written Rn = Rn,i + Rn,2, where 

M„ / j-l n \ M„ 

Rn,i=^an,j^iYl + - 2^^' i?n,2 = 2 a„,, [X , e] 

j=l \i=0 i=n-j+l j 1=1 



(51) 



1/2 

We shah show that Mn' R n converges in law, conditionally on X, to a normal distribution with 
variance 

i-l ry rl r-l 

War{e^) i i xh{x)h{y)dxdy + 8 {X, X) Var{e) / / h{x)h{y) min{x,y)dxdy, (52) 
Jo Jo Jo Jo 

and also that, conditionally on X, Rn/Mn"^ is asymptotically independent of {M^/n)^^'^(n in 
Corollary n in Section [l] Thus, in view of the results on the the pure noise and discretization terms 
in Corollaries n and 12 Theorem^ will then be shown. 

1/2 1/2 

To show this, we show in the following that M„ Rn.i and M„ Rn,2 are asymptotically normal 
given X, with mean zero and variances given by (|54() and ()57() . respectively. We then discuss the 
joint distribution of (M^/n)^/^Cn, Mn'^Rn,i and Mn'^Rn,2- 

Asymptotic normality of Rn^i- Once M„, < n/2, Write 

Af„-1 M„ Mn-1 Mn 

Rn,i = E 4 E ^ + E 4. E ^ - 2^^^- (53) 



Hence, 



j=0 j=i+l i=0 i=«+l 



M-1 M 



yar(MyXi) = 2M^Var{e')Y,iY. 

1=0 j=i+l 

= 2Varie'^) [ {[ h{y)dyf dx + o{l) 

Jo Jx 

/■I ry 

= AVar{e^) / / xh{x)h{y)dxdy + o{l) , (54) 
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while under Theorem [51 



E 



Mn -. i-1 n 



i=0 



2Ee2(l + o(M-^/2)) 



(55) 



1/2 

Since the Lindeberg condition is also obviously satisfied, I obtain that M„ Rn,i converges in law 
(conditionally on X) to a normal distribution with mean zero and variance given by equation (|54() . 

Asymptotic normality of the "cross term" Rn,2- As in the proof of Theorem |3J we proceed, 
without loss of generality, as if X were a martingale. Es in the proof of Theorem ^ we shall show 

1/2 

that Mn Rn,2 satisfies the conditions of t he version of the Mart ingale Central Limit Theorem 
which is stated in Corollary 3.1 (p. 58-59) of lHall and Hevdel h9H(k_ and calculate the as ymptotic 
variance. As in the earlier proof, we shall take, in the notation of iHall and Hevde ( 19801 ) . J-nj to 
be the smallest u-field making ej. , i = 1, and the whole Xt process, measurable. 



Note that, from ©, 



[X,e 



{n,K) 



n 



where 



-{Xt - 

(Xt . - Xt 



n.i — K ) 
nA~ K J 



{Xt - 



+K 



4=0 



X, 



if i = 0, • • • , - 1 

if i = K, - ■ ■ ,n — K 
if i = n-K + l,--- 



Thus, from (jSlf) . one obtains 



n M„, 



MyX2 = MV2j:e,5;^6(^]. 



(56) 



i=i j=i 



-1 /o 

Obviously, M„' Rn,2 is the end point of a zero mean martingale relat i ve to the filtration {J-n,j)- 
The conditional variance process (in Corollary 3.1 in iHall and Hevdf is given by (we use 

j A k = min(j, k)) 



n I M„ 



nE{e')E\E"-^'nl 



i=l \j=l 



nA n.i 



n Mn Mn 

MnVar{e)Y,Y.Y. , . 

1=1 j=i k=i 

n Mn Mn 

MnVar{e) ^^^.E '^'^{b'^f? + o,(l) 
i=i j=i k=i ■' 

Mn Mn 

0"n,j (l'n,k 



2MnVar{e)^Yl 



j=i k=i 



k 



(jAfc)[X,X](^-^^)+Op(l) 



1 fi 



JO 



h{x)h{y){x A y)dxdy {X, X) Var{e) + Op{l,, 



(57) 
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where remainder terms are taken care of as in the proof of Theorem |31 

By similar methods, the Lindeberg condit ion is satisfied ( c f. th e discussion in the proof of 
Theorem 11}. By Corollary 3.1 (p. 58-59) m iHaJI and HevHel (|l98nl ^ it follows that Mn^'^Rn is 



asymptotically normal (conditionally on X), with mean zero and variance given by (|57j) . This is 
what we needed to show. 

The joint distribution of {M'^/nY^'^C,n, Mn^'^Rn,! o-nd Mn^'^Rn,2- First of all, note that for all 



three quantities, we have satisfied the conditions of Corollary 3.1 (p. 58-59) of iHall and Hevde 



( 198nl ). This is with the exception of (their equation) (3.21), where we have instead used the 



Remarks following their corollary (and thus the convergence is conditional on X as opposed to 
stable with respect to the ci-field generated by both X and the etj- 

In terms of joint distribution, note first that the sum of conditional covariances (for each two of 
the three quantities {M^/ny/^Cn, Mn^'^R^,! and M^^i? 

n,2 converge to zero, by the same methods 
as above. In view of how Hall and Heyde's corollary implies their Theorem 3.2 (p. 58), the Cramer- 
Wold device now implies the joint normality of {M^/ny/'^(n, Mn^'^Rn,i and Mn^'^Rn,2, and also 
that they are asymptotically independent. Theorem 0] is then proved. 
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